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The effect of weak localization on spin relaxation in a two-dimensional system with a spin-split 
spectrum is considered. It is shown that the spin relaxation slows down due to the interference 
of electron waves moving along closed paths in opposite directions. As a result, the averaged 
electron spin decays at large times as 1/t. It is found that the spin dynamics can be described 
by a Boltzmann-type equation, in which the weak localization effects are taken into account as 
nonlocal-in-time corrections to the collision integral. The corrections are expressed via a spin- 
dependent return probability. The physical nature of the phenomenon is discussed and it is shown 
that the "nonbackscattering" contribution to the weak localization plays an essential role. It is also 
demonstrated that the magnetic field, both transversal and longitudinal, suppresses the power tail 
in the spin polarization. 
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Introduction 

The relaxation of non-equilibrium spin polarization is the central phenomenon in spin-dependent 
transport in semiconductor nanostructuresi One of the most efficient mechanisms of electron spin 
relaxation in III-V semiconductors is the well-known Dyakonov-Perel mechanism^ based on the clas- 
sical picture of angular diffusion of the spin vector in a random magnetic field. The field originates 
from the momentum-dependent spin-orbit splitting of the conduction band in the crystals with zinc- 
blend structured While passing through the crystal, the electron is scattered by impurities and its 
momentum changes randomly with time. As a consequence, the effective magnetic field also changes 
randomly with a correlation time of the order of the momentum relaxation time r. The spin relax- 
ation time ts is a characteristic time of the spin angular diffusion^ 1/ts ~ cf> 2 /r ~ £1 2 t , where 
4> ~ VLt 1 is the typical angle of the spin precession for the momentum relaxation time and f2 is 
the frequency of precession in the effective magnetic field proportional to the conduction band split- 
ting. For a two-dimensional (2D) case, when the electron motion in one direction is confined by the 
quantum well, the spin splitting and, hence, the precession frequency are proportional to the in-plane 
electron velocitj*^ Q, ~ v. So, at low temperatures, when inelastic processes can be disregarded, the 
spin relaxation rate for a 2D electron with a given energy E = mv 2 /2 is proportional to the particle 
diffusion coefficient 

— ~£>, D = v 2 t/2. (1) 

TS 

The effects of localization on the particle diffusion have been discussed in a great number of publi- 
cations. The first-order term in a series expansion of D in 1/fcZ (I — vt is the mean free path, k is 
the electron wave vector) is known as the weak localization correction^ ( for a review, see Ref. @), 
coming from the coherent enhancement of the backscattering amplitude. A remarkable feature of this 
correction is the logarithmic divergence at low temperatures in the 2D case. Equation Q implies a 
similar divergence of the spin relaxation rate. Such a divergence was first predicted by Singh£ in spin 
correlation functions for a system with spin-dependent impurity scattering. It was shown, however, 
that the quantum correction to the spin relaxation rate is not proportional to the quantum correction 
to the diffusion coefficient, as one might expect from Eq. (JTJ. A similar result was obtained in Ref. 
for a system with a spin-split spectrum. It was found in Ref. 10 that the weak localization slows 
down the spin relaxation of excitons in quantum wells, which leads to a 1/t power tail in the spin 
orientation. A similar effect was also predicted for electrons in 2D semiconductors with a zinc-blend 
crystal structure. 



In this paper, we consider the effects of localization on the spin relaxation for a 2D semiconductor 
with a spin-split spectrum. We show that the spin dynamics is described by a Boltzmann-type 
equation. In the first order in 1/kl, the localization effects can be taken into account by a nonlocal- 
in-time correction to the Boltzmann collision integral. This correction is expressed in terms of the 
spin-dependent return probability. We discuss the role of coherent returns at different scattering angles 
and show that the " nonbackscattering" contribution to the collision integral plays a key role. We solve 
the generalized kinetic equation and demonstrate that, at large times, the spin polarization decays as 
1/t. The magnetic field, both transversal and longitudinal, is found to suppress the long-living tail in 
the spin-relaxation. 



Derivation of the kinetic equation 

The Hamiltonian of a 2D with a spin-split spectrum is given by 

H = £- + ~u;cT + U(r), (2) 
2m 2 

Here p = pn is the in-plane electron momentum, in is the electron effective mass and <x is a vector 
consisting of Pauli matrices. The spin-orbit interaction is described by the second term, in which u> = 
u;(n) depends on the direction of the electron momentum u>i(n) = J2k n k^ki (i = x, y, z\ k — x, y). 
The matrix Cl = fiW + Cl^ is the sum of two terms: the so-called Bychkov-Rashba termi (with 

nonzero components Qxy = —Oyx ~ p) and &( 2 \ which is the Dresselhaus teiwm averaged over the 
electron motion along the z-direction perpendicular to the quantum well plane. The Bychkov-Rashba 
coupling depends on the asymmetry of the quantum well confining potential. Its strength can be tuned 
by varying the gate voltage^i The Dresselhaus term is present in semiconductors with no bulk inversion 
symmetry. The components of the matrix ti ^ are also linear in the in-plane electron momentum 

~ p and vary with well plane orientation with respect to the crystallographic axes^ (we neglect 
cubic Dresselhaus terms, assuming that the electron concentration is relatively small). We consider the 
scattering by the short-range impurity potential with the correlation function (U(r)U(r')) = jS(r—r'), 
where the coefficient 7 is related to the transport scattering time by r = h 3 /rwy. 

The classical spin dynamics is described by the kinetic equation^ For a homogeneous case, this 
equation is 

ds , . - ,„, 

— =w(n) xs + J s. 3 
at 

Here Jo is the Boltzmann collision integral and s = s(p, t) is the spin density in the momentum space, 
related to the averaged spin by S = J sd 2 p/(2irh) 2 . We assume that the spin splitting is relatively 
small: u>(n)T -C 1. This inequality provides r -C ts- The relationship between ts and the inelastic 
scattering time T- m varies with temperature. Here we focus on the case of low temperatures, assuming 
that ri n 3> ts- Then spins with different energies do not correlate with each other, and the solution of 
Eq. © at t > t yields S(t) = (m/2irfi 2 ) J SidE, where s ; = Si(E, t) = e^s^E, 0), Si(0, E) = (s(p, 0)) 
is the initial spin density averaged over the momentum direction, and T = 7\7 1 is the spin relaxation 
tensor (tensor of inverse relaxation times) given bj£ 

r lk = [s ik J2 Q si-J2 ^^ fc ] T l 2 - W 

8,1 I 

The conventional approach to the calculation of the correlations functions in weakly localized sys- 
tems is based on the Kubo formulaiSii An alternative approachi2iiiLii is to generalize the Boltzmann 
equation to include weak localization effects in the kinetic description. This approach may turn out 
to be more convenient when studying nonlinear and strongly nonequilibrium phenomena. To de- 
scribe quantitatively the weak localization phenomenon in the kinetic picture, one has to modify the 
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FIG. 1: Relevant irreducible diagrams a,b,c and the respective scattering processes a',b',c'. The Born 
collision process a is independent of the electron spin (its contribution is proportional to SapSe-f). Coherent 
backscattering b' ($ — $' » n ), as well as the processes c describing coherent scattering at an arbitrary angle 
(0 — 2n ), are spin-dependent due to rotation of the spin of an electron passing along the closed 

path. 

Boltzmann equation by introducing a nonlocal-in-time correction to the collision integra L 12 i 13 These 
corrections can be derived 13 from the diagrammatic structure of linear-response functions. Diagram- 
matically, the inclusion of a weak localization correction to the effective collision integral requires the 
consideration of the irreducible diagramsi&ii in Fig. ^y,c, in addition to the diagram for the Born 
scattering, shown in Fig. QJi. The crossed-ladder diagrams (lb) are usually considered to describe 
the coherent backscattering of the electron wave. A physical interpretation of the diagrams (lb,c) in 
terms of a small change in the effective differential cross-section for a single impurity was suggested in 
Ref.B It was based on an analysis of the interference contribution of trajectories propagating in the 
opposite directions along closed paths in terms of the phase stationarity requirement. This analysis 
shows that the diagrams (lb) correspond to a process shown in Fig.^/, which is indeed a coherent 
backscattering. Diagrams (lc) were found to describe coherent scattering processes with arbitrary 
scattering angles, shown in Fig.^J/. 

Next, we discuss the key points of a quantitative description of the weakly localized regime within 
the kinetic approach. We start with a brief discussion of the zero spin-orbit coupling. As can be seen 
from Figs.^'jC', the relevant processes contain the same closed paths, so the effective change in the 
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FIG. 2: Angular dependence of the effective cross-section modified by weak localization. The narrow peak 
at <£> — $' = 7r is due to the coherent backscattering shown in Fig. The enhancement of the coherent 

backscattering is accompanied by a reduction of scattering at other angles ( the process in Fig.^i'), the total 
cross-section being unchanged. 

differential cross-section of impurity 0, coming from both (lb') and (lc'), is expressed in terms of the 
return probability^*^ 



S 3M = M WU0)m-n)-l/2n] 

£>0 7T 



(5) 



Here So = 1/Nivt is the isotropic cross-section in the Drude approximation, Ni is the impurity 
concentration, A = 2ir/k is the electron wavelength, and W^(0) is given by 



W w (o) = ~ / ( //('-"n-(o.'i. 



(G) 



where W(0,t) = W(r, f)| r _>n is the probability density for a diffusing particle to return after the 
time t to the origin r = 0. The coefficient in the cross-section correction was found in Ref. 113 
by integration over small deviations of the electron trajectories from the trajectories (lb'), (lc') 
meeting the phase stationarity requirement. Physically, XI is the characteristic area of the region 
around the origin into which the diffusing electron should return for the effective interference to occur. 
The calculations show that the contributions of (lb') and (lc') have different signs. The positive 
contribution represented in Eq. (JSJ) by <5(<I> — 7r) comes from the process (lb'), while the negative 
one (the term — 1/2tt), from the process (lc'). In other words, the enhancement of the differential 
cross-section at the angle ir due to the coherent backscattering is accompanied by a reduction of the 
scattering in other directions, the total cross-section remaining unchanged (see Fig. |2J. We see that 
the correction to the effective impurity cross-section is w-dependent. Therefore, the correction to the 
collision integral in the time representation turns out to be nonlocal in time: 



5Jf(p, t) = (AZ/ttt 2 ) / dt' W(0, t-t') / d$'((S(* - $' - tt) - 1/2tt)/(p', t') 



(7) 



Here /(p, t) is the electron distribution function (we consider a homogeneous case, assuming that / is 
independent of r) and <&, $' are the angles of p and p'. Here and further we omit the "outflux" term 



-(AZ/ttt 2 ) f ^ dt' W(0, t-t')J d$'(<5($ - $' - tt) - l/2vr)/(p, t') in the collision integral. It makes 
a zero contribution, because the total cross-section is not changed by the weak localization. 

The probability density W u (r) is found as a sum over the paths involving different numbers of 
scattering events (see, for example, Ref. Ilfih 

W u (t) = J2wZ(t), (8) 



N 



where 



W2 



and 



(r) = J P(r - r N )P{r N - r w _i) . . . P(r 2 - r 1 )P(r 1 )dri . . . dr N , (9) 

P(t) = J— e -r/l+i"r/v / 1Q N 

At oj = 0, the function -P(r) is the classical probability density for an electron starting to move from r = 
to experience the first collision around the point r. In the framework of a diagrammatic approach, the 
function P(r) arises as the product of two spherical waves (retarded and advanced Green's functions) 
P(r) = jG% +hu G%. Here G R>A {r) = [Tif ,2 {m/h 2 )e ±lkr ' r ' 21 /V*^, k = y/^mE/h, and we took 
into account that e lkB + h " r e~ lkET w e lujr / v _ A path involving N scatterings contains (N + 1) functions 
Gr for the clockwise propagation along the path and (N + 1) functions Ga for the counterclockwise 
propagation. As a result, contains (N + 1) functions P(r). Using Eqs. ijHJl, (JHJ, JHJ, and l|10(l . one 
can show that in the diffusion approximation (lot C 1), when the typical number of scattering events 
along a path is large (N 3> 1), the function W(r, t) obeys the diffusion equation 

dW 

^- - DAW = 6(v)8(t). (11) 

Solving Eq. 111|) and taking r = 0, we find that the return probability is given by W(0,t) = l/4wDt, 
and W u (0) « (l/27r/ 2 )ln(l/rw). 

A generalization of the above results to a system with a spin-split spectrum is straightforward. 
Since the electron spin rotates while passing along a closed loop, electron Green's functions become 
operators with respect to the spin variables: 

™ ±ikr-r/2l 

G R ,A(r) = (T») 3/2 g r n , e~^^. (12) 

They are the products of spherical waves and spin-rotation matrices describing the electron spin 
precession. The precession frequency u?(n) depends on the propagation direction n = r/r. Due to 
the spin precession, the quantum correction to the effective cross section becomes spin-dependent: 
&S(<I>) — > 5SP~ <ae where the spin indices a, (3, 7, 9 correspond to the electron trajectories, as shown 
in Figs. Op', c': 

" - 1 ' = ~ W^ ae (Q)m - tt) - 1/27T], (13) 
The correction to the Boltzmann collision integral can be written as 

[SJffau)]^ = N iV [ SS^ ae (^ - *')/a«(p»d*', (14) 



where f a g(p' ,oj) is the momentum-dependent spin-density matrix. To derive the expression for W^ ia9 , 
we first introduce the probability density for a diffusing electron to arrive after N collisions at the 
point r with the spin rotated by an angle <fr 



(r, <t>) = 8{<t>- <f> N )P(r -r N )... P(r 2 - n)P(n)dn ...dr N . 



(15) 



The angle 4>n changes with the coordinates of the scattering points ri, . . . , rjy. One can find it from 
the matrix equation 

e -itj> N a-/2 _ e -ie(r-r N )o-/2 e -ie(r N -r 2 )o-/2 e ie(r2-ri)cr/2 e ~je(ri)cr/2 ^g-j 

where e(r) = o;(n)r/w. The spin-dependent return probability is then expressed via the total proba- 
bility density 

W u (r,ct>) = J2W?(r,<t>) (17) 

AT 

taken at r = 

W^ a8 {0) = J (/?|e-^ CT/ V)(0|e~ I0CT/2 |7}VMO,0) dA. (18) 

Here dA = g((j))d 3 (j> and function g(4>) is defined in Appendix 1X1 ( Note that, in the absence of the 
spin-orbit coupling, W U (Q, </>) = S(^))W u /g((()) and W^ a6 is expressed as W^ a6 = 5 a f}S 7 eW u ). What 
remains to be done is to find an equation for W u (r, cf>). To this end, the probabilities are related 
to each other by the recurrent equations 

5 (0)M^ +1 (r,0) = Jp( r ')W^(v-r',ct>')S(<l } -ct> , -A)dr'dA'. (19) 

The vector A = A r <0/ describes the change of the spin rotation angle in a ballistic path between two 
scattering. One can find it from the equation exp(—i(4>' + A)er/2) = exp(— ie(r')cr/2) exp(— i<//cr/2). 
Next, we sum Eq. (|19fl over N and take into account that A < 1. After combersome but straightfor- 
ward calculations, we find that W(r, <p,t) — t J W u (t, <fi) cxp(— iu}t)dio/2ir (the probability density to 
arrive after time t to the point r with the spin rotated by the angle cf>) is described by an equation 
similar to Eq. jP) : 

^~ D {j r -^f) W = 6(t)6(r)S(d>)/g(0). (20) 

Here L is the angular momentum operator acting on the functions of vector <p (see Ref. Il7h and 
(f2L) s ; = rjj^Lfc. The explicit expressions for L and for the common eigen- functions of L 2 and L z are 
presented in Appendix lAl Expanding <5(0) in a series over these functions, we keep the term with 
L = 0. The corresponding eigenfunction ^>o is independent of the angle </> and L^o = 0. This is the 
only term which survives at t 3> ts ■ The other terms decay exponentially with characteristic times of 
the order of t$- (This statement is not true for a degenerate case, when 17L depends on the component 
of L along a single axis. This case is discussed below). As a result, we obtain the following expression 
for the asymptotical behavior of 1^(0, <fi, t) 

W(0,d>,t) = -!—, for t»r s (21) 

Now, we write the distribution function as 

/ = If + sa, (22) 

where / is the particle density in the momentum space, related to the electron concentration by 
n = 2 J fdp/(2irh) 2 , I is the unit matrix, and s is the spin density. Substituting Eq. I|21l) and Eq. I|22f) 
into Eq. I|18|) and Eq. (|14|) , respectively, making a Fourier transform, and taking into account Eq. I|13l) , 
we obtain the weak-localization-induced correction to the collision integral 

(<5J/) = (XI/ttt 2 ) J dt'W(0,t-t') J d$'(«$($-$' -7r)-l/27r)/(p',0. ( 23 ) 
(XI/itt 2 ) f dt'W ik (0, t-t 1 ) f d$'(5(* - - tt) - l/27r)s fe (p', t'), (24) 



where Wik(0,t) and W(0,t) are given by 

Wi k (0,t) =J ^ lfe - 2e iefc sin 2 ^ W(0,&t)dA, W(0,i)= J coscf> W{0,4>,t)dA. (25) 
Here e = </>/</>. Using Eq. H21[). we find the asymptotical behavior of these functions 

W »^ = -£k> W ^ t) = -^Df fOT t>>TS - (26) 

Note also that the spin-orbit coupling can be neglected in the time interval r <C f <C . Hence, 
W(0, 0, i) ~ (5(0) and the expressions for Wik and W become 

WikM = 1^5i> W( °'* ) = ^rk' fOT T<<t «^ ( 2? ) 
We see that the difference between Eqs. (|2fi[) and Q27JI is in the numerical coefficients only 

Solution of the kinetic equation. The long-living tail in the spin polarization 

For the case with a spin polarization uniform in space, the generalized kinetic equation is 

Os S' — s 

— =u(n)xs + - \-8Js, (28) 

at t 

where 5 J s is given by Eq. I|24|l . This equation can be solved in the usual way£ Since w(n)r <C 1, the 
spin density can be represented as a sum of the isotropic part Si(£', t), which depends on the electron 
energy only, and a small anisotropic correction s a (p, t), which is linear in the electron momentum p : 

s = Si + s a . (29) 

Substituting Eq. <|29|l into Eq. Q28p. using of Eq. I|26() . and taking into account the equalities 

SJsi = 0, («yjs a ) = 0, (30) 

(here the angular brackets stand for averaging over the momentum direction) we obtain a closed 
relation for Sj 

§=-ff Si -j-r (3D 

m V tokij^ t-t'j y ' 

Assuming that the spin polarization was created at t = with a density Si(E, 0) and neglecting the 

quantum correction, we get the exponential relaxation Si(E,t) — 9(t)e~ rt Si(E,0) (here 8(t) is the 
theta- function) . This solution is valid until e -t / TS ~ 1/nkl. For larger times, the spin polarization 
should be found from the condition that the right-hand side of Eq. 131|) equals zero: S[(E,t) w 
(l/2nkl) fg dt'e~ Tt S[(E, 0)/(t — t'). So we find that the spin polarization has a long-living power tail 
at large times 

s i (E,t) = ^ M ^ rSi {E,0). (32) 

To conclude this section, we note that we neglected in our calculations the electron dephasing due 
to inelastic scattering. Such dephasing can be accounted for phenomenologically by introducing the 
factor exp(— t/r v ) into the right-hand side of Eq. I|32|) . Here t v is the phase-breaking time. 



The degenerate case 



Equation (|32|l is invalid for the degenerate case, when the spin precession frequency u>(n) is parallel 
to a certain vector u for any electron momentum: w(n) || u for any n. In the classical limit, the 
component of the electron spin parallel to u does not rela:»i£ T uu = and the two perpendicular 
components relax with equal rates Ti = T2 = T, the off-diagonal components of F being equal to zero 
(Ti v = = T±2 = ). This happens in symmetric quantum wells grown in the [110] direction, 5 
as well as in asymmetric quantum wells grown in the [001] direction, due to the interplay between 
Dresselhaus and Rashba couplings pliL 2 ^ To find the long-time asymptotic of the return probability, we 
write the formal solution of Eq. (|20[1 as 

W(Q,<t>,t)= I J^Le- D ^- ht '^ 2t 5{<t>)/g{0). (33) 



{2-Kf 

In the degenerate case, each of the three operators (OL) x , (£lL) y and (OL)^ is proportional to L u , 
which is the component of L along the precession axis. Therefore, these three operators commute 
with each other. Changing the integration variables q — > q iCt!L/v in Eq. (|33|) we obtain 

Thus, after travelling around a closed loop, the electron spin does not rotate at all. This can be 
interpreted as follows^ For the degenerate case, the spin rotation angle is simply given by <fi = 
J uxlt ~ J pdt. For a closed loop, we have J pdt = and, as a consequence, <p — 0. After substituting 
Eq. (|34|l into Eq. (|25|l and using Eqs. I|24[l. (|28() . one can see that the weak localization does not 
affect the longitudinal component of the spin density uSi(E, t) = usi(E, 0), while the relaxation of the 
perpendicular component is described by the following equation 

^--rf.-l/' Si ±u. (35) 



dt \ Tiki j_ oa t — t 

The integral in the right-hand side of Eq. I|35|l contains an additional factor 2 as compared with that 
in Eq. (|31|) . It can be shown that this factor arises from the contribution of the eigenfunctions with 
L = 1 to the long-time asymptotic (|34|l of the return probability. 21 The long-time asymptotic of the 
spin polarization is given by 

sdE,t) = — — si (.5,0), where s s (E,0) 1 u. (36) 
■kkI Tt 



Suppression of long-living polarization by the magnetic field 

Next, we consider the influence of the external magnetic field on the long-living power tail in the 
spin polarization. When an external field B is present, the spin rotation matrices appearing in Gr and 

G A become different: e -i">(«W2u e -i(«(n)+n )«rr/2t, for q r &nd e -iu(n)<rr/2v e -i(u)(n)-n )<Tr/2^ 

for Ga- Here, Flo is the frequency of the spin precession in the field B. So, Eq. H18|) is modified as 

w /*ra0 = J (/3| e -^/2| a )(0| e -^'-/2| 7 ) w u {Q,4>,(t>')dA dA', (37) 

where 

W u (r, 0, 0') = J2 1 5 i<t>- <t>N)5{<t>' - 4>' N )PB(r - r N ) . . . P B (r 2 - ri)P B (n)dri . . . dr N . (38) 

N J 



Her. 
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P B (r- r') = P(r-r')exp x r ']j , (39) 

and is the component of the magnetic field normal to the quantum well plane. The vectors 0at and 
<p' N should be found from Eq. (|16[l with e(r) = [u>(n) + flo]r/v for 0jv and e'(r) = [<*>(n) — fio]f*A> 
for cj)' N . Writing out the recurrent equations for the consecutive terms in Eq. I|38[) . we can derive the 
following expression 

^ - iO (L L')W - D (^-i^- » " (L + L,) ) W = 5(t)5(r)6(4>)5(d>')/9 2 (0h (40) 

for the probability density W(r,<j>,<j>',t) — r J W w (r, 0, 0') exp(— iu)t)dw/2TT. Here the operators L 
and £/ are given by Eq. IjAlfl (with the replacement ~ * 4>' f° r L')- The term —i2eA/ch in Eq. (|4U[) 
is responsible for the magnetic field effect on the electron orbital motion. 22 Here A = [Bj_ x r]/2 is 
the vector potential. We assume that the magnetic field is small 

MoTs < 1, and DeB ± T S /hc < 1. (41) 

For B = 0, the solution of Eq. gUJ) is 

W(r, 0, 0', *) = W (r, 0, t) 5(0 - <f>')/g((t>), (42) 

where W(r, <f>, t) obeys Eq. H20JI . Using Eq. I|42|l. we rewrite the long-living solution H21(l (which 
corresponds to zero total angular momentum J = L + L') in terms of two angles 0, 0': 

-^ 6{ *J ] = for 5 = 0. (43) 

AnrDt g(</>) 7 (2tt) 2 ff (</») V j 

Here 7(g) = Dq 2 is the eigenvalue of the operator 7 = Z?(q — ClJ /v) 2 at J = and 5{<p — cf>')/g(<p) is 
the respective eigenfunction. For B ^ 0, the long time dynamics of the spin relaxation is determined 
by the eigenvalues of the operator 

7= -ifl Q {L - L') + D (-id/ dr - 2eA/ ch- ClJ . (44) 

First, we neglect the term — ifio(L — L'). In this approximation, the eigenvalues of 7, corresponding 
to J = 0, are given by 

7» = 7J.(»+l/2), (45) 

where 7j_ = 4eB±D /He and n = 0, 1, 2, Since the operator D(—id/dr — 2eA/ch) 2 has a discrete 

spectrum, we have to make the replacement 

^_ e -7(,)* ^ V - ^ (46) 

(2tt) 2 47r£> „ 47T J Dsinri(7 X t/2) v ; 

in Eq. I|43[l . The second step is to take into account the term —iCto(L — I/), taking it to be a 
small perturbation. In the first order of the perturbation theory, this term leads to the mixing of the 
eigenfunctions with J — 1 to the eigenfunction with J — 0. This mixing can be disregarded at f^o <C T. 
Corrections to the eigenvalues 7„ arise in the second order only. They are calculated in Appendix [H] 
Using Eq. I|B5|) . we can show that the spin polarization dynamics at large times is described as 



Thus, we see that the magnetic field does suppress the long-living tail in the spin polarization. 



Discussion 



Next we discuss briefly the physical meaning of the results obtained. Our calculations were based on 
the interpretation of the weak localization phenomena in terms of two scattering processes: coherent 
backscattering (see Fig.^ 3 ') and coherent scattering at an arbitrary angle (see Fig.H]c'). The existence 
of the long-living spin polarization can be explained as follows. Both coherent scattering processes 
were shown to be nonlocal in time. In other words, the transition of a spin s to a spin s', caused 
by coherent scattering, takes a certain time t, which may be relatively long: t 3> t$. The coherent 
scattering events do not change the direction of the electron spin. Indeed, as seen from Eq. 1(26(1 . 
Wik(0, t) ~ 5ik and, as a consequence, s' || s. Therefore, the electrons involved in such a scattering 
can keep memory about the initial spin polarization during the time much longer than rg. The power 
law 1/t is due to the proportionality of the probability of coherent scattering to the probability of 
diffusive return. 

The collision integral accounting for both coherent processes does not change the total scattering 
cross-section. Therefore, 

5 J fi = 0, 6 J Si = 0. (48) 

Using Eqs. 123(1 . 124(1 . 126(1 . and taking into account that / a (— p, w) = — / a (p, w) and s a (— p,u>) — 
— s a (p, u>), we obtain 

6 J ,Up,u) = / a (p,w), 6 J s a (p,w) = - ln 2 ( ^ r T) s a (p, W ), for ur s « 1. (49) 

Eqs. ((48(1 . 1 (49(1 imply that the effect of localization on the angular spin diffusion, as well as that on the 
particle diffusion, can be accounted for by the w-dependent renormalization of the transport scattering 
time. However, the quantum corrections to this time have differnt signs for the particle and angular 
spin diffusion: 



1 _ l -ln(l/cr)/2^ 



Ttr 



for the particle diffusion and 



1 + In (1/W) /2irkl 

(51) 



for the spin diffusion. This implies that the Eq. Q is invalid in the quantum case, or, more precisely, 
it relates ts and r t ' r rather than t$ and r tr . 

An important comment should be made concerning the role of the coherent nonbackscattering 
contribution (processes d in Fig. Neglecting this contribution, we have SJsi ^ 0. It can be 
easily shown that this leads to a physically meaningless result for the spin relaxation rate. Therefore, 
the correct treatment of the effect of weak localization on the spin relaxation is only possible when 
the nonbackscattering coherent processes are taken into account (the role of such effects for particle 
diffusion was discussed in Ref. fl4f) . It is worth noting that the weak localization effects are usually 
considered to be due to the coherent backscattering only. The point is that the quantum correction 
to the conductivity is usually calculated by means of the Kubo formula, which expresses conductivity 
in terms of the current-current correlation function. This approach focuses on the calculation of the 
velocity correlation function, which depends on the anisotropic part of the distribution function / a ; 
so, there is no need to know corrections to The situation is quite different for spin relaxation which 
is due to the relaxation of the isotropic part of the distribution function^ 

To conclude, we have discussed the long-time dynamics of the spin polarization in a 2D disordered 
semiconductor. It is shown that, at large times, the spin relaxation slows down due to weak localization 
effects. An analytical expression for the long-living tail of the spin-polarization has been derived. 
The magnetic field, both transversal and longitudinal, suppresses this tail, restoring the exponential 
relaxation. 
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APPENDIX A: 



The analytical expression for the angular momentum operator L is 



i ^ ^ 

1 cot — 

2 2 



e x 



e x 



d_ 
~d4> 



2 X 90 



(Al) 



with e = 4>/<fi. The components of this operator obey the usual commutation rules, [Li, Lj] = ieijiLi. 
The common eigenfunctions of L 2 and L z are Wigner's rotation matrices^ Dj( 1M ,(4>): 

L 2 D L MM ,{4>) = L(L + l)D L MM ,(4>), L z D L MM ,{ct>) = MD L MM ,{cj>). (A2) 

The values of L may be integer and half-integer L = 0, 1/2, 1, 3/2, . . . (as usual, M = —L, . . . , L). 
There is also (2L + 1) degeneracy with respect to M'. The eigenfunction corresponding to L = is 
equal to unity 

* = D° 0Q = 1, L* = 0. (A3) 

The orthogonality conditions arc 17 

Sll^mm^m'M'/^L + I). (A4) 

Here the integration is made over all possible transformations of the SU (2) group. The integration 
measure is taken in the invariant form: 

dA = 9 {<t>)d 5 4>, where g(<f>) = ■ (A5) 

Note that SU(2) transformations are usually parameterized by Euler angles < a < 2ir, < (3 < 
7r, — 2ir < 7 < 2ir . For such a parametrization, the invariant integration measure is given byi£ 
dA = sin/3 da d(i d^/lQ-K 2 (the expression for L via the Euler angles is also given in Ref. Hit) . It 
is convenient for us to use the components of vector <p (0 < <\> < 2ir) instead of Euler angles. This 
gives 24 dA = g((j))d 3 <p. The expansion of the delta function 5(<fi — </>') in Wigner's functions is 

5(0-0')= E (2^ + 1)^W(0)^mV(0'M0)- (A6) 

L,M,W 

The solution to Eq. H2(J[) can be represented as a sum over all angular momenta L = 0, 1/2, 1, 3/2, .... 
However, as follows from Eq. (|18fl only two terms (L = and L = 1) contribute to the spin-dependent 
return probability. Indeed, the action of the operator L on the spin-rotation matrices is given by 

L(/?|e-"^/ 2 |a) = ^Ml(/?'| e -»W 2 | a ), ( A7 ) 
L(P\e-^ 2 \a)(6\e~^ 2 \ 7 ) = ^'^+<W<t^ (p\ e -»~/*\a){0'\e- i *'/ 2 \<y). (A8) 

As follows from Eq. (|A8(I . the projection of the operator L onto the subspace formed by products of 
two rotation matrices is given by 

L , = +* (Ag) 

where erW and cr^ are the Pauli matrices acting on the first and second rotation matrices, respec- 
tively. Therefore, the angular momentum L' can only be or 1 (singlet and triplet contributions). 

Using Eq. J20J) and the property J dA* 1 L^ , 2 = - J dA^ 2 i'^i, valid for arbitrary functions ^i(4>) 
and ^2(0)) we can see that the function 

W^ a6 {r,t) = J </3|e-^ CT/2 |o ; >(6'|e-^ CT / 2 | T )V^Cr, 0, i) dA (A10) 

obeys the equation similar to Eq. I|2U|) 



d Id iHL' 
dt \ dr v 



W 0' ja e' = S (t)S(r)5 aP S ie , (All) 



P/3'88' 

where matrix elements l^'na/ggi are given by Eq. I|A8(1 . Note that Eq. (|Alip was used in Refs. 

for calculation of weak localization corrections to conductivity. The alternative derivation of the 

operator L' was given in Ref. [27] 



APPENDIX B: 



Using Eq. I|A6|) . we expand Eq. (gSJ as 

5{<t> - 0') 



L Mi=-L 



where 



(Bl) 



*f lM2 (0,0') = V2TTT £ ^ Mi (0)^mm 2 (0') 



(B2) 



M=-L 



is the full set of the eigenfunctions (L = 0, 1/2, 1, 3/2, . . . ; Mi = — i, . . . , L; M% = —L, . . . ,L) for 
zero total angular momentum j-qj^ IlM2 = 0. In the second order of the perturbation theory, the term 
— zf2o(L — L') leads to the corrections to the eigenvalues of the operator 7. These corrections are 
different for the functions \j/^ flMl with different values of L. They are calculated as 



D -l ( ^AfiMi 



flo(L-L') (foj/v) n (L-L' 



(B3) 



Next, we diagonalize the operator (flj/v) 2 in the subspace formed by the three functions (L 2 
L' x )^f lM '\ (L y - L' y )^ M \ and [L x - L' Z )^ M \ Direct calculation gives 



A 7i = ^^W-no. 



The expression for the long-living solution becomes 

1 7± /2 



W(0,<f>,<f>',t) 



AttD sinh(7 ± t/2) 



^2 V2L~+l e~ AlLt y, * 



Mi Mi 
L ' 



(B4) 



(B5) 



where L = Q, 1/2, 1,3/2, .... Substituting Eq. (jB5jl into Eq. ^3 and using Eqs. O, O, and 
we derive Eq. (|I7ll . 



